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We obtain sufficient conditions for the efficient simulation of a continuous variable quantum algo- 
rithm or process on a classical computer. The resulting theorem is an extension of the Gottesman- 
Knill theorem to continuous variable quantum information. For a collection of harmonic oscillators, 
any quantum process that begins with unentangled Gaussian states, performs only transformations 
generated by Hamiltonians that are quadratic in the canonical operators, and involves only measure- 
ments of canonical operators (including finite losses) and suitable operations conditioned on these 
measurements can be simulated efficiently on a classical computer. 
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Quantum mechanics allows for information processing 
that could not be performed classically. In particular, 
it may be possible to perform an algorithm efficiently 
on a quantum computer that cannot be performed effi- 
ciently on a classical one. Significant effort is now un- 
derway to construct quantum algorithms and processes 
that yield such a speedup. The Gottesman-Knill (GK) 
theorem Q for discrete-variable (qubit) quantum infor- 
mation provides a valuable tool for assessing the classical 
complexity of a given process. Essentially, it states that 
any quantum algorithm that initiates in the computa- 
tional basis and employs only a restricted class of gates 
(Hadamard, phase, CNOT, and Pauli gates), along with 
projective measurements in the computational basis, can 
be efficiently simulated on a classical computer. (For a 
precise formulation and proof of this remarkable theo- 
rem, see |Q, page 464.) The GK theorem reveals that 
a large class of quantum algorithms do not provide a 
speedup over classical processes. In fact, recent work has 
placed even stronger constraints on the potential speedup 
of fermionic quantum computers || . 

In addition to the successes of qubit-based algorithms, 
quantum information over continuous variables (CV) 
has yielded many exciting advances, both theoretically 
and experimentally, in fields such as quantum tclepor- 
tation j|, quantum cryptography |@, 0, H [|, and 



potentially quantum computation 1 10 . CV algorithms 



could also perform computational tasks more efficiently 
than is possible classically. To assess the computational 
complexity of these tasks, it is necessary to develop an 
extension of the GK theorem: what continuous variable 
processes can be efficiently simulated on a classical com- 
puter? As a CV quantum information process involves 
coupled canonical systems, this question of efficient clas- 
sical simulation is related to asking under what condi- 
tions a quantum mechanical system can be modeled by a 
classical one. As noted by Feynman 111}] , a key advantage 



of a quantum computer is its ability to simulate quantum 
systems that cannot be efficiently simulated classically. 

The issue of efficient classical simulation of a CV pro- 
cess is more involved than for the discrete case. One 
notable difference is that the quantum states and the 
unitary transformations involved are described by real- 
valued (as opposed to integer-valued) parameters, and 
these parameters must be described on a discrete classical 
computer with some assumption of error or limited pre- 
cision. Also, the states used in CV experiments are ap- 
proximations to the idealized computational basis. These 
basis states are infinitely squeezed states whereas any ex- 
perimental implementation will involve finitely squeezed 
states thus deviating from their idealized form [lPf, A 
good classical simulation must be robust against such de- 
viations. Measurements are part of the quantum compu- 
tation and, even in the computational basis, are subject 
to experimental constraints (such as photodetection effi- 
ciency). Classical simulation must also incorporate these 
measurements. 

Despite these complications, we prove in the follow- 
ing an extension of the GK theorem for continuous vari- 
ables; i.e., we present a set of sufficient conditions for 
a CV quantum information process which, if satisfied, 
ensure that it can be efficiently simulated on a classical 
computer. To prove this theorem, we employ the tech- 
niques of stabilizers [^) that are used for qubits. Using 
the stabilizer formalism, it is often possible to simulate a 
quantum information process by following the evolution 
of a set of operators, the Pauli operators, rather than the 
evolution of quantum states. For CV processes, we show 
that it is more natural to analyze stabilizers in terms 
of the algebras {i.e., Hamiltonians) that generate them, 
rather than the groups themselves. We define analogs of 
the Pauli and Clifford algebras and groups for CV and 
construct sets of gates (as unitary transformations) that 
can efficiently simulate any arbitrary transformation in 
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these groups. Any algorithm or process constructed out 
of these Clifford group transformations can be efficiently 
modeled by following the evolution of the Pauli operators 
rather than the states of the system. 

The standard Pauli group Q n for CV quantum compu- 
tation on n coupled oscillator systems is the Heisenberg- 
Weyl group HW(rt), which consists of phase-space dis- 
placement operators for the n oscillators. Unlike the dis- 
crete Pauli group for qubits, the group HW(n) is a con- 
tinuous (Lie) group, and can therefore only be generated 
by a set of continuously-parameterized operators. The al- 
gebra hw(n) that generates this group is spanned by the 
2n canonical operators qi, pi, i = 1, . . . , n, along with the 
identity operator /, satisfying the commutation relations 
fej^j] = ihSijJ. For a single oscillator, the n = 1 alge- 
bra is spanned by the canonical operators {q,p, 1} which 
generate the single oscillator Pauli operators 
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with q,p G R. The Pauli operator X(q) is a position- 
translation operator (translating by an amount q), 
whereas Z(p) is a momentum boost operator (kicking 
the momentum by an amount p). These operators are 
non-commutative and obey the identity 



X(q)Z(p) = e-^Z(p)X(q) . 



(2) 



On the computational basis of position eigenstates 
{|s); sel} fLOL pi p|, tlle Pauli operators act as 



X(q)\s) = \s + q), Z{p)\s) 



exp(-ps)\s) 



(3) 



Note that it is conventional to use highly squeezed states 
to approximate position eigenstates; these states satisfy 
the orthogonality relation (s\s') = 5(s — s') in the limit 
of infinite squeezing. 

The Pauli operators for one system can be used to 
construct a set of Pauli operators {Xi(qi), Zi(pi); i = 
1, . . . , n} for n systems (where each operator labeled by 
i acts as the identity on all other systems j ^ i). This 
set generates the Pauli group Q n . Note that the Pauli 
group is only a subgroup of all possible unitary transfor- 
mations. It is not possible to construct an arbitrary uni- 
tary transformation using only the Pauli operators X(q) 
and Z(p); the Pauli group only describes transformations 
generated by Hamiltonians that are linear in the canon- 
ical variables. 

For issues of classical simulation, we will be interested 
in transformations that lie in the Clifford group. The 
Clifford group N(Q n ) is the group of transformations, 
acting by conjugation, that preserves the Pauli group Q n ; 
i.e., it is the normalizer of the Pauli group in the (infinite- 
dimensional) group of all unitary transformations. 

Theorem 1: The Clifford group N{Q n ) for con- 
tinuous variables is the semidirect product group 



[HW(n)]Sp(2n, R), consisting of all phase- space transla- 
tions along with all one-mode and two-mode squeezing 
transformations. This group is generated by inhomoge- 
neous quadratic polynomials in the canonical operators. 

Proof: The most straightforward method to iden- 
tify the Clifford group will be to identify its algebra. 
The Clifford algebra consists of all Hamiltonian oper- 
ators H c satisfying [Hh w ,H c ] G hw(n) for all Hi lw <E 
hw(n). This algebra must obviously include the alge- 
bra hw(n), and thus hw(n) is a subalgebra of the Clif- 
ford algebra. In addition, this algebra includes all ho- 
mogeneous quadratic polynomials in the canonical oper- 
ators {cji, Pi', i — 1, . . . ,n}. This algebra of quadratics 
consists of Hamiltonians that generate one-mode squeez- 
ing transformations [for example, the Hamiltonian Hs = 
5 (qP + P<i)] i an d also interaction Hamiltonians that gen- 
erate two-mode squeezing transformations (for example, 
the interaction Hamiltonian H mt = qi <S>p2). The alge- 
bra of homogeneous quadratic polynomials in the canon- 
ical operators is known as the linear symplectic algebra 
sp(2n,K). 

Together, the algebras hw(n) and sp(2?i, M) form a 
larger algebra, consisting of inhomogeneous quadratic 
Hamiltonians in the canonical operators {qi,pf,i — 
1,... ,n}. This algebra is the semidirect sum algebra 
[hw(n)]sp(2n, R), with hw(n) as an ideal. The group gen- 
erated by this algebra is the semidirect product group 
[HW(n)]Sp(2n,M). This group includes phase-space dis- 
placements (the Pauli group), as well as the squeezing 
transformations (both single- and two-mode) of quantum 
optics @. (QED) 

In order to describe a quantum information process 
as a circuit, it is necessary to find a set of transforma- 
tions (gates) that generate the Clifford group; these gates 
will serve as building blocks for arbitrary Clifford group 
transformations. Following the derivation by Gottesman 
et al flj|, a set of gates will be defined in terms of the 
elements of the Clifford algebra (i.e., the Hamiltonians) 
that generate the transformations. 

The SUM gate is the CV analog of the CNOT gate 
and provides the basic interaction gate for two oscillator 
systems 1 and 2; it is defined as 



SUM = exp(--^i <&p 2 



(4) 



This gate is an interaction gate operation on the Pauli 
group Q2 for two systems. Referring to the definition (Q) 
for the Pauli operators for a single system, the action of 
this gate on the Q2 Pauli operators is given by 



SUM : Xx(q)®I 2 

h®X 2 {q) 
h®Z 2 {p) 



X 1 {q)®X 2 {q), 
Z x (p)®h, 
h®X 2 {q), 
Z x (p)- X ® Z 2 {p) . 



(5) 
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This gate describes the unitary transformation used in 
a back-action evasion or quantum nondemolition pro- 
cess Q. 

The Fourier transform F is the CV analog of the 
Hadamard transformation. It is defined as 

F = exp(lJ(<f +p 2 )), (6) 

and the action on the Pauli operators is 

F : X(q) -» Z(q) , 

Z{p)^X{p)-\ (7) 

The 'phase gate' P(rj) is a squeezing operation for CV, 
defined by 

P(V) = exp(^ 7?( 7 2 ) , (8) 

and the action on the Pauli operators is 

P( V ) :X(q)^e^ 2 X(q)Z(r ] q), 

Z{p) - Z{p) . (9) 

(The operator P{rj) is called the phase gate, in analogy 
to the discrete- variable phase gate P |H| , because of its 
similar action on the Pauli operators.) 

For discrete variables, it is possible to generate the 
Clifford group using only the SUM, F, and P gates [ p"5| . 
However, for the CV definitions above, the operators 
SUM, F, and P(r?) are all elements of Sp(2n,R); they 
are generated by homogeneous quadratic Hamiltonians 
only. Thus, they are in a subgroup of the Clifford group. 
In order to generate the entire Clifford group, one re- 
quires a continuous HW(1) transformation [i.e., a linear 
Hamiltonian, that generates a one-parameter subgroup 
of HW(1)] such as the Pauli operator Xiq). This set 
{SUM, F, P(t}), X(q);r),q e K} generates the Clifford 
group. 

We now have the necessary components to prove the 
main theorem of this paper regarding efficient classical 
simulation of a CV process. We employ the stabilizer 
formalism used for discrete variables and follow the evo- 
lution of the Pauli operators rather than the states. To 
start with, let us consider the ideal case of a system with 
an initial state in the computational basis of the form 
|<7i, <?2, ■ ■ • j ?n)- This state may be fully characterized 
by the eigenvalues of the generators of n Pauli opera- 
tors {qi, fa, ■■ ■ , q n }- Any continuous variable process 
or algorithm that is expressed in terms of Clifford group 
transformations can then be modeled by following the 
evolution of the generators of these n Pauli operators, 
rather than by following the evolution of the states in 
the Hilbert space C 2 (M. n ). The Clifford group maps lin- 
ear combinations of Pauli operator generators to linear 
combinations of Pauli operator generators (each qi and pi 
is mapped to sums of qj , pj , j = 1, . . . , n in the Heisen- 
berg picture). For each of the n generators describing the 



initial state, one must keep track of 2n real coefficients 
describing this linear combination. To simulate such a 
system, then, requires following the evolution of 2n 2 real 
numbers. 

In the simplest case, measurements (in the computa- 
tional basis) are performed at the end of the computa- 
tion. An efficient classical simulation involves simulat- 
ing the statistics of linear combinations of Pauli operator 
generators. In terms of the Heisenberg evolution, the qj 
are described by their initial eigenvalues, and the pj in 
the sum by a uniform random number. This prescription 
reproduces the statistics of all multi-mode correlations 
for measurements of these operators. 

Measurement in the computational basis plus feed- 
forward during the computation may also be easily simu- 
lated for a sufficiently restricted class of feed-forward op- 
erations; in particular, operations corresponding to feed- 
forward displacement (not rotation or squeezing, though 
this restriction will be dropped below) by an amount pro- 
portional to the measurement result. Such feed-forward 
operations may be simulated by the Hamiltonian that 
generates the SUM gate with measurement in the com- 
putational basis delayed until the end of the computation. 
In other words, feed-forward from measurement can be 
treated by employing conditional unitary operations with 
delayed measurement ^], thus reducing feed- forward to 
the case already treated. 

In practice, infinitely squeezed input states are not 
available. Instead, the initial states will be of the form 

Si(n)(g)5a(r2)<8)---(g)5 n (r n )|0 ) 0,... ,0), (10) 

where |0) is a vacuum state and S{r), r 6 K is the 
squeezing operation which can be expressed directly in 
terms of elements of the Clifford group. Now the vacuum 
states may also be described by stabilizers {qi+ifti, 92 + 
ip2, ... , q n + ipn} which are complex linear combina- 
tions of the generators. Combining the initial squeezing 
operators into the computation, a classical simulation re- 
quires following the evolution of 4n 2 numbers (twice that 
of infinitely squeezed inputs due to the real and imagi- 
nary parts). Measurements in the computational basis 
are again easily simulated in terms of this Heisenberg 
evolution, by treating each of the qi and pi as random 
numbers independently sampled from a Gaussian distri- 
bution with widths described by the vacuum state. Sim- 
ulation of measurement plus feed-forward follows exactly 
the same prescription as before. 

The condition for ideal measurements can be relaxed. 
Finite efficiency detection can be modeled by a linear loss 
mechanism jl6| . Such a mechanism may be described by 
quadratic Hamiltonians and hence simulated by the Clif- 
ford group. Note that the Clifford group transformations 
are precisely those that preserve Gaussian states; i.e., 
they transform Gaussians to Gaussians; this observation 
allows us to remove our earlier restriction on feed-forward 
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gates and allow for classical feed-forward of any Clifford 
group operation. Note that non-Gaussian components to 
the states cannot be modeled in this manner. 

Finally, it should be noted that modeling the evolution 
requires operations on real- valued (continuous) variables, 
and thus must be discretized when the simulation is done 
on a discrete (as opposed to analog) classical computer. 
The discretization assumes a finite error, which will be 
bounded by the smaller of the initial squeezing or the 
final detector 'resolution' due to finite efficiency, and this 
error must remain bounded throughout the simulation. 
As only the operations of addition and multiplication are 
required, the discretization error can be kept bounded 
with a polynomial cost to efficiency. 

Thus, we have proved the extension of the GK theorem 
for continuous variables: 

Theorem 2 (Efficient Classical Simulation): Any 
continuous variable quantum information process that 
initiates with Gaussian states (products of squeezed dis- 
placed vacuum states) and performs only (i) linear phase- 
space displacements (given by the Pauli group), (ii) 
squeezing transformations on a single oscillator sys- 
tem, (Hi) SUM gates, (iv) measurements in position- 
or momentum- eigenstate basis (measurements of Pauli 
group operators) with finite losses, and (v) Clifford 
group [HW(n)]Sp(2n, R) operations conditioned on clas- 
sical numbers or measurements of Pauli operators ( clas- 
sical feed-forward) , can be efficiently simulated using a 
classical computer. 

We could summarize the conditions (i-iii) by simply 
stating (i-iii) transformations generated by Hamiltonians 
that are inhomogeneous quadratics in the canonical oper- 
ators {fa, pi, i — f , . . . , n}, which is equivalent. Thus, 
any circuit built up of components described by one- 
or two-mode quadratic Hamiltonians [such as the set 
of gates SUM, F, P(t]), and X(q)}, that initiates with 
finitely squeezed states and involves only measurements 
of canonical variables may be efficiently classically simu- 
lated. 

As with the discrete-variable case, these conditions 
do not mean that entanglement between the n oscilla- 
tor systems is not allowed; for example, starting with 
(separable) position eigenstates, the Fourier transform 
gate combined with the SUM gate can lead to entan- 
glement. Thus, algorithms that produce entanglement 
between systems may still satisfy the conditions of the 
theorem and hence may be simulated efficiently on a 
classical computer; included are those used for CV quan- 
tum teleportation Q, quantum cryptography |7|, s|. p |, 
and error correction for CV quantum computing |12, [13 1. 
Although these processes are of a fundamentally quan- 
tum nature and involve entanglement between systems, 
this theorem demonstrates that they do not provide any 
speedup over a classical process. Thus, our theorem pro- 
vides a valuable tool in assessing the classical complexity 
of simulating these quantum processes. 



As shown in JlO| ] , in order to generate all unitary trans- 
formations given by an arbitrary polynomial Hamilto- 
nian (as is necessary to perform universal CV quantum 
computation), one must include a gate described by a 
Hamiltonian other than an inhomogeneous quadratic in 
the canonical operators, such as a cubic or higher-order 
polynomial. Transformations generated by these Hamil- 
tonians do not preserve the Pauli group, and thus cannot 
be described by the stabilizer formalism. Moreover, any 
such Hamiltonian is sufficient [ [To| ]. One example would 
be to include an optical Kerr nonlinearity |l7| , but there 
is a lack of sufficiently strong nonlinear materials with 
low absorption. Alternatively, it has recently been pro- 
posed that a measurement-induced nonlinearity (using 
ideal photodetection) could be used in an optical scheme 
without the need for nonlinear materials in the compu- 
tation |]l5[ The physical realization of such nonlin- 
carities is an important quest for quantum information 
theory over continuous variables. These nonlinear trans- 
formations can be used in CV algorithms that do not 
satisfy the criteria of this theorem, and which may pro- 
vide a significant speedup over any classical process. 
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